RANK-CRANK TYPE PDES AND 
GENERALIZED LAMBERT SERIES IDENTITIES 



SONG HENG CHAN, ATUL DIXIT AND FRANK G. GARVAN 
Dedicated to our friends, Mourad Ismail and Dennis Stanton 

Abstract. We show how Rank-Crank type PDEs for higher order Appell functions 
due to Zwegers may be obtained from a generalized Lambert series identity due to 
the first author. Special cases are the Rank-Crank PDE due to Atkin and the third 
author and a PDE for a level 5 Appell function also found by the third author. These 
two special PDEs are related to generalized Lambert series identities due to Watson, 
and Jackson respectively. The first author's Lambert series identities are common 
generalizations. We also show how Atkin and Swinnerton-Dyer's proof using elliptic 
functions can be extended to prove these generalized Lambert series identities. 



1. Introduction 

F. J. Dyson [9], [TOj p. 52] denned the rank of a partition as the largest part minus 
the number of parts. Dyson conjectured that the residue of the rank mod 5 divides 
the partitions of 5n + 4 into 5 equal clases thereby providing a combinatorial inter- 
pretation of Ramanujan's famous partition congruences p(5n + 4) = 0(mod5). He 
also conjectured that the rank mod 7 likewise gives Ramanujan's partition congru- 
ence p(7n + 5) = (mod 7). Dyson's rank conjectures were proved by A. O. L. Atkin 
and H. P. F. Swinnerton-Dyer [3]. The following was the crucial identity that Atkin 
and Swinnerton-Dyer needed for the proof of the Dyson rank conjectures. It was first 
proved by G.N. Watson 
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Throughout we use the standard g-notation 
(x) := (x;q) := 1, 

n— 1 

W.. ' . 

m=0 

j ■Em\n • [^1 j " ' ' j Xmi q\n • Q/ -^1) 5 -^m; -£m> ?)n- 



(a?;9)n := JJ(1 



when n is a nonnegative integer. Assuming \q\ < 1 we also use this notation when 
n = oo by interpreting its meaning as the limit as n — > oo. Later M. Jackson 
proved an analogue of the above identity, 

C 2 [C 2 ]oo[<]ook/C]oo y (-l)V"("+i)/g | [CUC 2 UxCUC/xUq)l 



C, l^Joolt, Joo , 1 x„ 5n(n+l)/2 / , X ) Q q^) 



x / /^-5n /-5n+5 \ 

Ef_lf>(»+l)/2 / __S _J> \ 
1 ; y \l-<-V 



Recently, the first author [HI p. 603] found a generalization of the above two identities, 
namely, 

X?[X 2 /X 1 , ■■■ , X m jx U X X X m , ■■■ , X ± X 2 , Xftao ^ (_ 1 )n g (2m+l)n(n+l)/2 



[^i]oo[x 2 , • • • , Xm} 2 ^ 1 - 

Xi/x 2 , ■ ■ ■ , Xi/x m , Xi, XiX m , , X\X 2 , X^]oo 

00 



\Z / X\ : Z J X 2: • • • 1 zj X m , Z, ZX m , , ZX\ 



oo 



X 2 [X2/S3, • • • ,£2/£ 



-(2m+l)n (2m+l)(n+l) ' 



x {-iTq^^ 2 7 w/ + ^ + idem(x 2 ; x 3 , 

n z~^o \ 1 — zq n /x 2 1 — 2;x 2 g / 

00 

£ (-I)' 



00 / -(2m+l)n (2m+l)(n+l) ' 
\n (2m+l)n(n+l)/2 ( J^l , ^1 



1 — zq n jx\ 1 — zx\q n 
where g(a\, a 2 , ■ ■ ■ , a m )+ idem(ai; a 2 , ■ ■ ■ , a n ) denotes the sum 

n 

(?(aj, a 2 , • • • , cti-i, di, a>i+i, ■ • • , a. m ), 

in which the i-th term of the sum is obtained from the first by interchanging a\ and a^. 

Equation (11. 3p was proved using partial fractions. Indeed, the m = 1 case of (11. 3p 
is equivalent to (11.11) . while the m = 2 case is equivalent to (II. 2p . The fact that the 
right-hand side of (11.21) is independent of x, and that the right-hand side of (11.31) is 
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independent of x 2l x 3 , ■ ■ ■ ,x m seems to be intriguing at first. Indeed, one purpose of 
this article is to show that the left-hand sides of (II. 2p and (11.31) are really elliptic 
functions of order less than 2, in fact entire functions as we show, in the respective 
variables (x for ( 11.21) and x 2 for ( 11.31) while holding £3, • • • ,x m fixed) and therefore 
that they must be constants which are nothing but the right-hand sides of (II. 2p and 
(II. 3p respectively. Since ( II. 2p follows from (II. 3p . we show this only for (II. 3p . This is 
done in Section [2j 

Let N(m, n) denote the number of partitions of n with rank m. Then the rank 
generating function R(z, q) is given by 



00 n 



^^yjyjjvK^^yj^-j^. (1 . 4) 

In PP, G. E. Andrews and the third author defined the crank of a partition, a partition 
statistic hypothesized by Dyson in [9]. It is the largest part if the partition contains 
no ones, and otherwise is the number of parts larger than the number of ones minus 
the number of ones. For n > 1, we let M(m,n) denote the number of partitions of n 
with crank m. If we amend the definition of M(m,n) for n — 1, then the generating 
function can be given as an infinite product. Accordingly, we assume 

M(0, 1) = -1, M(-l, 1) = M(l, 1) = 1, and M(m, 1) = otherwise. 

Then the crank generating function C(z, q) is given by 

C( „)^5; % .)A-^. (1.5) 

n=0 m=—n x ' v ' 

Atkin and the third author [2] found the so-called Rank-Crank PDE, a partial differen- 
tial equation (PDE) which relates R(z, q) and C(z, q). To state this PDE in its original 
form, we first define the differential operators 

6, = z^, 6 q = q±. (1.6) 
Then the Rank-Crank PDE can be written as 

z(q)UC%z,q)] 3 = U. + ^ + ^sA R*(z,q), (1.7) 



where 



R*(z,q) :-- 
CT(z,q) :-- 



R(z,q) 
l-z ' 

C(z,q) 
l-z 



In [2], it was shown how the Rank-Crank PDE and certain results for the derivatives 
of Eisenstein series lead to exact relations between rank and crank moments. As in 
define Nk(m,n) by 

^N k {m,n)q n = -— ^(_l)n-l g n((2fc-l)n-l)/2+|m|n( 1 _ ^ (1.9) 
n>0 n=1 
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of any positive integer k. When k — 1 this is the generating function for the crank, and 
when k = 2 it is the generating function for the rank. When k > 2, N k (m,n) can be 
interpreted combinatorially as the number of partitions of n into k — 1 successive Durfee 
squares with fc-rank equal to m. See p3J Eq.(l.ll)] for a definition of the A;-rank. We 
define 

n 

R k (z,q) :=J2Y1 N k (™,n)z m q n . (1.10) 

n>0 m=—n 

From [T3], Eq.(4.5)], this generating function can be written as 

^nf+n|H hn|_j 

Rk ^ q)= £ 77) rT7(o) (zQ ) ' ( L11 ) 

when k > 2. In Section [31 we show that R k (z,q) is related to the level Ik — 1 Appell 
function 

tf'-'i, 2>():=s (li — . (1 . 12) 

n=— oo 

We obtain the following 
Theorem 1.1. For k > 1, 

Rk(z,q) 
1 



k-3 



Jfe-1, 



(?)c 

where 



1 - zJE^- 1 )^, 9 ) - ^i, 2fc -i(g) + z(l -z)J2 z m Q 2M ^{q) 



m=0 



(1.13) 

e^-M = E (-i)v ((2fc ~ 1)n+j)/2 , (i.i4) 



/or j = 1, 3, . . ., 2k -3. 

This theorem generalizes Lemma 7.9 in [12] which gives a relation between the rank 
generating function R(z, q) and a level 3 Appell function. The k = 1 case of the theorem 
gives the familiar partial fraction expansion for the reciprocal of Jacobi's theta product 

(z)oo(z- l q)oc m P- 1], H3 P- 136], since E"=-oo(- 1 ) n 9" (n+1)/2 = °- 

A few years ago the third author found a 4 th order PDE, which is an analogue of the 
Rank-Crank PDE and is related to the 3-rank [13]. To state this PDE we define 

G( 5 )( M ):~S( 5 '( M ). (1.15) 

Then 

24(q)UC*(z,q)} 5 
= 24(l-10$ 3 (g))G^(z,g) 
+ (1005, + 505* + 1006 q 5 z + 355 2 z + 205 q 5 2 z + 1005* + 10<5, 3 + 5f) G (5 \z, q), (1.16) 
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where 

00 3 n 

•»(«) = E£^=- < L17 > 

n=l 

This PDE can be written more compactly as 

24(g) 2 0O [C*(z, q)f = (Hi - E A ) G^(z, q), (1.18) 
where H* is the operator 

H* := 5 + 105, + 55 z + 5 2 z , 

and 

E 4 :=E 4 (q) := 1 + 240$ 3 (g), 

is the usual Eisenstein series of weight 4. The PDE (11.161) was first conjectured by the 
third author and then subsequently proved and generalized by Zwegers [21] . It was also 
Zwegers who first observed that (11. 1 6 [) could be written in a more compact form. We 
now describe Zwegers's generalization. Define for I G Z >0 , the level I Appell function 

as 

Mu, v) := Mu, v; r) := z 1 ' 2 £ (L19) 

z — ' 1 — zq n 

n=— 00 

where z = e 2nm , w = e 2mv , q = e 2mT , and define the modified rank and crank generating 
functions as follows. 

1/2 —1/24 

K := K(u;t) := R(z,q), (1-20) 

1 — z 

1/2 —1/24 

C := C(u; r) := - ^— C(z, q). (1.21) 

Here and throughout we assume Imr > so that \q\ < 1. Then the following theorem 
due to Zwegers gives for odd I, the (I — l) th order analogue of the Rank-Crank PDE. 

Theorem 1.2 (Zwegers|21j). Let I > 3 be an odd integer. Define 

Id Id 2 l(2k-l) 

(2m) 2 du 2 12 2 ' 

:= ^2fc-i"H2fc-3 • • •'Hz'Hi, (1-22) 

where E 2 (t) = 1 — 24 Y^=i o~i(n)q n is the usual Eisenstein series of weight 2 with 
o~a{ n ) — ^2d\n^ a - Then there exist holomorphic modular forms fj (which can be con- 
structed explicitly), with j = 4, 6, 8, • ■ ■ ,1 — 1, on SX 2 (Z) of weight j , such that 

(«-5)/2 \ 

K [a-i)/ 2 ] + f^jtw A( u ,o) = (i- ly.rfe 1 , (1.23) 

k=0 J 

where r] is the Dedekind r]-function, given by r)(r) = q 1 ^ 2i (q)oo- 
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Zwegers proved this theorem using the formulas and methods motivated from the 
theory of Jacobi forms. In contrast to this, the proof of the Rank-Crank PDE by Atkin 
and the third author, which corresponds to the / = 3 case of Zwegers's PDE, depends 
upon simply taking the second derivative with respect to ( of both sides of (II. ip . The 
main goal of this paper is to show how a generalized Rank-Crank PDE of any odd order 
follows from the Lambert series identity (II. 3p in a similar fashion. We will obtain 
Zwegers's result in a different form. In our form the coefficients are quasimodular 
forms rather than holomorphic modular forms, but in contrast, our coefficients are 
given recursively. See Theorem 14.41 and Corollary 14.51 

This paper is organized as follows. In Section [2J we prove (II. 3p using the theory 
of elliptic functions. Then in Section [3] we prove Theorem II. 1| which is the theorem 
that relates Rk(z,q) with the level 2k — 1 Appell function £( 2fc_1 )(,2, q). In Section H] 
we prove our main result that shows how (ll.3p can be used to derive the higher order 
Rank-Crank-type PDEs of Zwegers. 

In the light of (11.191) . it should be observed that the identities (11.11) and (II. 2p are really 
the identities involving certain combinations of level 3 and level 5 Appell functions 
respectively while (11.31) is an identity involving a combination of level (2m + 1) Appell 
functions. However, the analogue for level 1 Appell function which cannot be derived 
from (11.31) was found by R. Lewis fTJ, Equation 11] and is as follows. 

[^]oo[C ]oo(<?)oo \ r / \n ra(n+l)/2 

[zCUCUzC-'U ) q Vi - zrK i - <<t 




2. General Lambert series identity through elliptic function theory 

Atkin and Swinnerton-Dyer's proof of (II. ip depends in essence on the theory of 
elliptic functions. In this section, we show how this method of proof can be used to 
prove (II .3p . Let 

z = e 2 ™, xi = e 2 ™\ x 2 = e 2mw , (2.1) 

and let 

X] =e 2ma \ j = 3,---,m, (2.2) 

where u, v, w, a^, . . . , a m are all complex numbers. Also recall that q = e 2nlT , where Im 
r > 0. 
Let 

[6] oo =: J(a, q) =: J(a), where b = e 2nia , a E C. (2.3) 

Then using the Jacobi triple product identity [U p. 10, Theorem 1.3.3], we easily find 
that, 

„-„7ria „-l/8 

J(a,g) = ^ 0(a), (2.4) 



where 



6(z) = 6(z;t) := e K«+lA+ 2 K"+i) H) . (2.5) 



n=— oo 
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Comparing this with the classical definition of 9i(z) [Tlj p. 355, Section 13.19, Equation 
10], we find that upon replacing q by q 1 ^ 2 in this classical definition, 9(z) = —9 1 (z). 
From [2QI P- 8] 5 we see that 

9(z + l) = -9(z), (2.6) 

e( z + t) = -e- mr ~ 2mz 6{z), (2.7) 

9(-z) = -9(z), (2.8) 

#'(0;r) = -27^/8^. (2.9) 

Using (12.61) and (12. 7p . we have 

J{a + l,q) = J{a,q), (2.10) 

J(a + r,q) = -e- 27Tia J(a,q), (2.11) 

J(a -T,q) = -q~ 1 e 2,Ria J{a, q), (2.12) 

J(-a, q) = - e - 27Tia J(a, q). (2.13) 

Using (12.31) . we rephrase (11.31) as follows: 



f™™J(w - V)J(W + V) J(2v) yr Jjoj - v) Jjoj + v) ^ (_ 1 )n g (2m+l)n(n+l)/2 



J(v — w) J(v) J(v + w) J(2v)(q) 2 Kl t-i- J{v — dj)J{v + dj) 
J(u — v)J{u — w)J(u)J(u + w)J(u + v) 3< ~.^ m J(u — a,j)J(u + dj) 

, f 2m(v-w) J(v)J(2v) yr J (v - Oj ) J (v + a, ) 

1 J{w)J{2w) 3 ll m J(w - a,) J( w + a,) 

X V (_l)n g (2m+l)n(n+l)/2 / _ + £ + idem( W ; a 3 , • • • 

n=— oo ^ ' 

00 / -2vm;(2m+l)n 2™(2m+l)(n+l) \ 

E(_n« (2m+l)n(n+l)/2 ( _£ , t ) (2 14) 

\ / "* \ 1 — e 2 7ri(u-i))Qn ^ — e 2ni( u + v )q n ) ' 

n=— oo ^ ' 

Fix 03, • ■ ■ , a m , consider the left-hand side of (12. 14ft as a function of w only and denote 
it by g(w). Let fi(w) denote the expression in line 1 of (12. 141) . f^iw) the expression in 
line 2 of (12. 14ft and /s(w) the expression in lines 3 and 4 of (12. 14ft . Then, using ( I2.10p . 
( l2TTj) and (|2TT2|) . we see 

/i(w + 1) = /iW = /i(w + t), 
/ 2 («J + l) = / 2 H = / 2 (w + r), 

/ 3 (w + l)=/ 3 H. (2.15) 
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Another application of (12. lip and (I2.12p gives 
fs(w + r) 

_ e 27Ti( - v - w -^J{v)J{2v) -j-r J{v - aj )J(v + dj) 

~ J{w + r)J{2w + 2r) 3 1 | m J(w + r - aj )J(w + r + a,) 

00 / „-27ri(w+r)(2m+l)n „27ri(u>+T)(2m+l)(n+l) \ 

x C_l)n (2m+l)n(n+l)/2 / £ j_ £ \ 

/ j \ / " I g27ri(n— ui— t)qU J g27rj(u+«'+T) ( jn i 

n=— 00 ^ y ■■ ^ / 

, y^ c 2^-a fc ) ^(^)J(2^)J(^-^-r)J(t; + w + r) -j-r J(t; - flj -) + a 3 -) 
£^ J(a k )J(2a k )J(a k -w- r)J(a k + w + r) 2< ^ +1 J(afc - aj)J(a k + aj) 

j^k 



00 / -27rja fc (2m+l)n 2izia k (2m+l)(n+l) \ 

x C_l)n „(2m+l)n(n+l)/2 / ± y £ \ 

Z-^i V 7 y V.1 — e 2iTi(u-a k )qn J _ e 2iri(u+a k ) qn J 



X 



n=— 00 

27rCT+47rimw J(v)J(2v) TT J(f — aj)J(v + 0j) ( /_-.\n (2m+l)(n+l)(n+2)/2 

JH J(2 W ) 3 ll m J( W - % ) J(«, + a,) V n £^ J 

„-27rilt)(2m+l)(n+l) „-(2m+l)(n+l) 00 „27iW2m+l)n „(2m+l)n \ 
« 1 y^ (2m+l)n(n-l)/2f V \ 



~y £,2Tri(u—w) qfi 



! y^ e 27ri( P -a fc ) J(^)J(2i;)J(^-w)J(^ + w) -r-r 7(^-^)7^ + ^) 
^ J(a k )J(2a k )J(a k - w)J(a k + w) 2< fJ^ +1 J(a k - a j )J(a k + aj) 

3^k 



00 / -2nia k (2m+l)n 2izia k (2m+l)(n+l) \ 

x \^ (_l)n (2m+l)n(n+l)/2 ( , £ ] 

V / y \1 — g27ri(u-a fc )gn J _ e 2ni(u+a k ) q-a J 



n=— 00 



= /sH- (2.16) 

Thus from (12.1 5p and (I2.16p . we deduce that g is a doubly periodic function in w with 
periods 1 and r. Our next task is to show that g is an entire function of w and hence 
a constant (with respect to w). We show that the poles of g at w = u and w = —u are 
actually removable singularities by proving that Xmi w ^ t ± u {w =F u) (f2{w) + h{w)) = 
which readily implies that lim w _>± u (w =F v )g{w) = 0. Let 

A := A(v,a 3 , ■ ■ ■ ,a m ;q) := J(v)J(2v) JJ J(u - a,) J(u + a,). (2.17) 

3<jr'<m 

Next, applying |OJ, (EZE5D and (JH, we see that 
lim (w - u) (f 2 (w) + f 3 (w)) 

= A\im(w-u)\- Av-^Av + rvMlo rr 1 

w^u y J(u — w)J(u + w)J(u — v)J{u)J{u + v) 3 ^.^ m J{u — aj)J{u + aj) 

e 2wi(v-w) ^ / J 

3<j'<m 



_i_ TT [ t i_ (_-\\n (am+l)n(n+l)/2 

JH J(2w) 11 J(w - aj)J(w + a,) \1 - e 2 ^ u ~ w ) ^ 1 ; y 



n— — oo 
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-2mw(2m+l)n 00 2niw(2m+l)(n+l) ' 



00 27ri«;(2ni+l)(n+l) \ ^ 

i (_iy» „(2m+l)n(n+l)/2f \ I 

T 1 I q I _ e 27ri(u+w)o»» i J 



n=— 00 



1 / -^ 1/8 (g)L J_ 

J(lt)J(2u) J(U- 0^)7(14 + Oj-) V ^'(0) 2?TZ 

= 0. (2.18) 

Similarly, \im w ^_ u (w + u) (/^(w) + fs(w)) = 0. Now the only other possibility of a pole 
of g is at 0, which arises from fi and f 3 each having a pole at 0. Again, to show that 
this is a removable singularity, it suffices to show that \im w ^ w (fi(w) + fs(w)) = 0. 
To show this, we need Jacobi's duplication formula for theta functions [TSJ p. 488, Ex. 
5] 

e(2w)e 2 e 3 e 4 = 29(w)e 2 (w)e 3 (w)9 4 (w). (2.19) 



Let 



d r>( n . \ ._ r 2Kimv J ( 2v ) TT J(a j -v)J(a j + v) (~l) w g 

B - B(u,v,a 3 ,--- ,a m ,q) .-e — — [[ — — ^ TZ 



n (2m+l)n(n+l)/2 



J(v) 11 (J(aj)) 2 ^ l_ e 2niu q n 

3<j<m v \ 311 n=-oo * 

(2.20) 

Then from (gJl and fl2T2D|) . 

lim w (A(w) + h(w)) 

w— >0 

f J(w -v)J(w + v) . e 2ni ^- w ^ -,-r 1 
lim w < B — - — rr-r h A— — — - 

,,,-4.(1 I ( T T n,A 1 9in\ 11 



w^o [ (J(w)) 2 J(w)J(2w) ^ m J(w-a 3 )J(w + aj ) 



00 / c -2iriw(2m+l)ri 2-Kiw(2m+l)(n+l) \ ~\ 

x (_l)n (2m+l)n(n+l)/2 / _^ , £ \ I 

/ j \ / " y j g27ri('U— w)qn 'Y £,2ni(u+w) qn J J 



n=— 00 

m 



I lim^ Vc 2 ^- afc) ^)^)^-^)^ + ^) TT J(^-aj)J(^ + Qj) 
w-M> ^ J(a k )J(2a k )J(a k -w)J(a h + w) 2< fJ^ +i J(a k -a j )J(a k + a j ) 

w f 9{w-v)9{w + v) 2m 1/4f , e- 5m n 1 
«-o*(«;)\ ' 9 ^e(2 W ) 3 ll m J( W -a J )J( W + a J ) 

00 / -27riw(2m+l)n 27Uio(2m+l)(n+l) \ "\ 

E/_ 1 xn (2m+l)n(n+l)/2 / _£ , f ) I 

V / " ^ g27ri(M— w)qn ^ Q2ni(u+w) qn J J 



X 

n=— 00 



1 >to#^, (2-21) 



0'(O) ^0 0(w) 
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where 

D(w) := D(w, v, a 3 , ■ ■ ■ , a m ; q) := B6(w - v)6{w + v)- Ae 2mv q^^f^w) 

n/ , N e- 5niw 6(w) -,-r 1 

E(w) := E[w; u, a 3 ■ ■ ■ , a m ; q) := — { [ 



^(2^) 35 ^ m J(w - %)J(w + aj) 



00 / „-27rilu(2m+l)ri 27mu(2m+l)(n+l) \ 
x (_l)n (2m+l)n(n+l)/2 / _i_ , £ ] 



n=— 00 



(2.22) 

Now using fl2~19D . f[2T4l and (12351) . we find that as w -> 0, 

D , X_,_„ fl3r X e 2^ Mg l/4 (g) 2 o « ( _ 1) n g (2 m+ l)n(n + l)/2 

lU ' J W (-l)^2 e -3brf(as+".+« ra )(J( ag )... J( am ))2jL^ !_ e 2™ g n 

= 0, (2.23) 

which is observed by putting the expressions for B and C back in the first expression 
on the right side in f |2.23[) . Thus, 

\im w {h{w) + h{w)) = ^\. (2.24) 
w-x) v v ' v " 6»'(0) 

Now we calculate -D'(O). 

D'(w) = B (e'{w - v)6(w + v) + B{w - v)9'(w + vj\ - e 2 ™ Atf^qf^Ef (w). (2.25) 
Using (J23D and (gISD , we have 

7riui(2m+l) /3 /) /) i 
%) = (-l)-2 g ¥(g)V- 2)i W TT ' 



20 2 H0 3 («;)04H 3 < A ^ m 0(™ - a,)^^ + a,) 



00 / -2vri«;(2m+l)n 27riuj(2m+l)(n+l) \ 

x V (_i)» g (am+i)K«+i)/2 _ + £ . (2.26) 

/ -< y q2-ki(u—w) qn J q2tti(u+w) q-n J x ' 



n=— 00 



Differentiating both sides with respect to w and simplifying, we obtain 

£'H = i(-l) m - 2 g I ^(g) 2 ( m - 2 )e-™( 2m+1 ^ 2 03^4 

•/o m I i\ I flgM g^(w) / e'jw-aj) 9\w + a 3 ) 

1 + j + e 2 (u,) + 3 («O *s(«0 3 <^ m V 0(™ - %) 0(™ + «,) 



X 



3<j'<m 

x / -2niw(2m+l)n 2niw(2m+l)(n+l) \ 

E/m (2m+l)n(n+l)/2 / _£ , f \ 
V / " q2-ki(u—w) qn J Q2iri(u+w) qn J 



n=— 00 



+ 2 H0 3 M04W 3 |J m - o,)0(u; + a,) F H ^ (2 ' 27) 
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where 

00 / p -2niw(2m+l)n 2niw(2m+l)(n+l) \ 

F(w) := F(w, u, m; q) := V 2 m +i)n(n+i)/2 + £ 

n=— 00 ^ ' 

(2.28) 

It is straightforward to see that 

00 f_1^n „(2m+l)n(n+l)/2 

F'(0) = 2K( 2m + 1) ^ [ V-e^a" " (2 " 29) 

n=— 00 

From (12. 8p . we have 

e'(-z) = 9'{z). (2.30) 

Then letting w ->• in (E2ZD, and using (EHD, (ESQ]), ( 12T29D and the fact that 0^(0) = 
for 2 < k < 4, we find that 

£'(0) = 0. (2.31) 

Using fl2T30|) and f[23T]) in fl2T25l) . we finally deduce that L>'(0) = 0. 

With the help of (I2.24p . this then implies that lim w ^ w (f\(w) + fziw)) = and 
thus lim^-i.o wg(w) = 0. Thus w — is also a removable singularity, which implies that 
g is an doubly periodic entire function and hence a constant, say K (which may very 
well depend on v). Finally, since J(0) = 0, we have 

00 / e -2niv(2m+l)n e 2™(2m+l)(n+l) \ 

K = g(v) = 2J ("I) <? ( + ) ( + ^ ( — e27ri(u _„ ) + ^ _ e27r i( u +t,) n J ' 

71= — CO ^ ' 

This completes the proof. 

3. Proof of Theorem 11.11 
From [131 Eq.(4.3)], we see that 

i 00 / 1 — 1 n \ 

= -^^- i r- l v n{{2k - l)n - l)/2 (i - ?») (t ± ^ + T^^T 



n=l 
i 00 



1 — « r 

( n((2fc-l)n+l)/2 1 g 



(<Z)oo„t^ 1-^-V 

(3.1) 



Replacing 2; by 2; 1 in (II. lip and (13. ip . we see that 
Rk(z,q) 



Z / 1 \n-l„n((2fc-l)n+l)/2 ^ 9™ 



n((2fc-l)n+l)/2 



Z /_-|\n-l n((2fe-l)n+l)/2 / 1 _ (1 ~ z)g n \ 



(<?) 



n— — 00 
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/ 



(?)c 



(?)c 



E (-irv 

Ti — — OO 

Z' oo 

1 - E c- 1 )"?" 



((2fc-l)ra+l)/2 



(i---) E 



__-|\»lg7i((2fc-l)n+3)/2 



n — — oo 



1 — zq n 



((2fe-l)n+l)/2 



+ E 



-l) n q 



n n((2fc-l)n+3)/2 



V 



n— — oo 



1 — zg" 



-Z01,2Jfc-l(?) ^ 



(9)00 (g)c 
-^i, 2 fc-i(g) . 2 



(00 

n#0 



-l) n g 



n n((2fc-l)n+3)/2 



1 — zg r ' 



(?)c 



(<?)c 



1 E(-i)v 



((2fc-l)n+3)/2 



n — — 00 
n#0 



fc— 2q(k— 2)n j (zq n ^~ 



1 - zg* 



+ 



1 — zg r 



-z0i )2 fc-i(<z) ^ 



(9)00 (g)c 



-^i,2fc-i(g) + ^ 



i + (i-,) E(-!) 



n^n((2fc-l)n+3)/2 



n — — 00 
n#0 



r fc-2„(fc-2)n fe - 3 

zq 1 



t^— + E 



m ran 
£ q 



m=0 



(?)c 



(g)c 



1 + 2 



fe-2 



E 



n (2fc-l)n(n+l)/2 



-l) n g 



n— — 00 

A;— 3 00 



m=0 n=-oo 



-Z0l,2k-l(q) z 



(9)00 (g) c 



1 - zg n 

+ (1 _ 2 ) ^ z m ^ (_l)"g™(( 2fc - 1 )™+ 2m + 3 )/ 2 
m=0 

00 

71= — OO 

3 00 \ 



v / _^«g(2A:-l)n(n+l)/2 



m=0 n=— 00 



/ fc— 3 

^- ( -^i, 2fc -i(g) + 2 fe - x (l - g) + ^(1 - ^) E 



This completes the proof of Theorem 11.11 



2m+3,2fc-l 



(?) • 



m=0 



4. Higher order Rank-Crank- type PDEs 



In this section we show how the generalized Lambert series identity ( 11. 3D can be used 
to derive general Rank-Crank PDEs of the type found by Zwegers. 

4.1. The idea. First we let = £*, 1 < i < m in (11. 3p to obtain 

m— 1 

*UC, *, ?) (?) 2 oo = W(C, *, ?)+E ?) 5Wf.i(C i+l , *, g)-^o |W (C, 9) £ (2m+1) (^ g), 

3=1 

(4.1) 
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where 



00 / A-fcn fk(n+l) \ 

S fe (C, z, q) := £ (-l)"/" ( " +1)/2 i _\ + , (4.2) 

n=~oc \ ' 

for k odd and 

F ,UC,q):=C [ ^r^, (4-3) 



j 00 



|>-(m-l) /•- (m-2) . . . /— (m-j)l 

^(c^)- 1 — [c : +2; . . . ; cm+ ; +1L J °° (fon<j<m-i), (4.4) 

[a— (m— 1) A— (m— 2) a— 2 1 /■ f-1 Am ATO+l] 

' ^ ' [<" m , ^C"^- 1 ), • ■ ■ , <" 2 , Z(~\ Z, Z(, Z(\ ■ ■ • , zC m ~\ Z^L ' ^ ' 1 

The basic idea is to apply the operator D 2m to both sides of (14. ip where 



5 W 



We will also need the differential operator 

%l := k5 z + 2k5 q + 5 2 z . (4.7) 

We note that the operator 7-L* k differs from Zwegers's Tik although they are similar. First 
we need to write the functions Sw(z, q) and Sk(z, q) as double series. Throughout we 
assume that < \q\ < 1, z {q n : n e Z} U {0} and C £ {z ±l q n : n G Z} U {0}. We 
obtain 

OO / OO OO 

5 fc (C, q) = J^(-i) V n(n+1)/2 C kn Yl z m c m q mn + C k(n+1) z m Cq r ' 

n=0 \ m=0 m=0 

OO / OO OO 

^ V -Qn^fcn(n— 1)/2 j ^fcra ^ ^ ^—m^rn^nn _|_ ^fe(— n+1) ^ ^ £~ m (^" m g r ' 

n=l \ m=l m=l 

OO OO 

^ ^ ^ ^ ( ^w^m^fen(n+l)/2+mn ^—kn—m _|_ £fc(n-t-l)+m^ 

n=0 m=0 

OO OO 

^ ^ ^ ^ ( ]^ n 2 _ Tnqkn(n— l)/2+mn ^Afen+m _|_ a— kn+k— ^ 

and 

£^(2, g) = ^(-l) n z ra /"' n+1 " 2+ ™ — ^(-l)VV'"" 1 ' , ' 2+mn . (4-9) 



n=l m=l 



OO OO OO OO 



n=0 m=0 n=l m=l 



We have 

Theorem 4.1. Suppose k is odd and 1 < £ < k — 1. T/ien 

DtS k (C,z,q) = P kji (Ht)EW(z,q), (4.10) 



14 SONG HENG CHAN, ATUL DIXIT AND FRANK G. GARVAN 

where 

^H£f^=^- (4.1D 

Proof. Suppose k is odd and 1 < i < k — 1 . First we prove that 



P k ,l{x) = {\k - \ Vk 2 + Axf + {\k + \ Vk 2 + Axf. (4.12) 

We have 

, , / p\ 

{\k - ±Vk 2 + 4x) e + (|ife + | Vk 2 + Axf = 53 I 2 • ) k l ~ 2j 2 l ~ l (k 2 + 4x) J 

i=o ^ J ' 

= 53 53 ( 2 •) ( J )x m ^- 2m 2 2m - m = 53 53 ( .) r ) x m ^- 2m 2 2m - m . 

j=0 m=0 V^/ V m / m=0 I J=m \2j/ \ m J J 

The result (I4.12p now follows from the binomial coefficient identity 

^\2})\m (<-2m)!m!' 



which we leave as an exercise. 

We observe that if x = km + m 2 + k 2 n(n + 1) + Imnk then 

A; 2 + Ax = (k + 2m + 2A;n) 2 , 

|A; — ^Vk 2 + 4x = — /cn — m, 



+ i\/fc 2 + 4x = + 1) + m, 

and we see that 

Similarly we find that 

^ ^fcn+m + ^-fcn+fc- m ^ = (jfc n+m )^ + (_jb n+ jb_ m )< = P ke (-km+m 2 +k 2 n{n~l)+2mnk) . 
We note that 

Ul ( q ^(n+l)/2+mn z m^ = ^ + m 2 + + ^ + 2mnA; ) ^(n+lj^m^mj ^ 

Thus 

n /„tn(ti+l)/2+tnn jm (/ — kn—m , /-kn+k+m\\ d lu*\ ( „kn{n+l) / 2+mn „m\ 

L>e\q K " z {(, +(, ))=r k/ (H k ){q z ), 

and 

n / „kn(n+l) / 2+mn „m r /-kn+m , /-—kn+k—m\\ r> /„b(n-l)/2+mn„-ra\ 

Utyl z \s )) = PkA H k) [Q z )■ 

The result f HTLOl) follows from equations (@~8D and (l4~9]h □ 

Next we calculate D 2m of each term in (14. ip . 
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4.2. The term Y m ((, z, q). It is clear that the function Y m ((, z, q) has a zero of order 
2m at ( = 1. It is well-known that 

2m 

£> 2m (/(C)) = ^5(2m,j)/ 0) (l), 
i=i 

where the numbers S(2m, j) are Stirling numbers of the second kind. Since S(2m, 2m) = 
1 it follows that 

D 2m (Y m ((,z,q)))=Yi 2m \l,z,q) = (-l)™" 1 ^)! (m+1)! (m-1)! [C*(z, q)} 2m+1 (q) 2 ^ 1 

(4.14) 

by an easy calculation. 

4.3. The term F 0tm ((,q). By logarithmic differentiation we have 

S C Fo,m(C,q) = L , m ((,q)F , m ((,q). (4-15) 

where 

L , m (C, g) = K ,m(C) - (m + 1) E ^_ Cm+lgi " Y^^lJ 

OO x 

+ m E(i 

1=1 x 



- C"V 1 - C~"V 

= K , m (0 -(m+i) (c +1 qr - cr™ - vr + ™ E (cvr - (cvr 

i,n>l i,n>l 

(4.16) 

^o,m(C) = m + J m (C) - J m _i(C), (4.17) 

m 

E< n 

■/m(C) = ^ • (4-18) 

n=0 

For any positive integer k we define 

i 00 „2fc-l n r> 

G 2fc := G 2k (q) := -((1 - 2fc) + £ = + $ 2 *-i(<z), (4.19) 

n=i 

where B 2n is the (2n)-th Bernoulli number, and 

OO 

$2*-i := $2fe-i(?) := E ff 2fe-i(n)g n . (4.20) 

n=l 

The function G 2 fc is a normalized Eisenstein series. For k > 1 it is an entire modular 
form of weight 2k. We need the following 



Lemma 4.2. If m and a are positive integers, then 

hi 
+ 



D a (U0) = f^y ((m + l) a+1 - 1) . 
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Proof. Suppose m and a are positive integers. It is well-known that 

e x -l ^ k\ y ' 

k=0 

By taking the logarithmic derivative of (£ m+1 — l)/(£ — 1) we find that 

JmiC) = m + (m + I)- — . 

Hence by (14.2 li we have 

Bi 



k=0 

The result now follows since 



J m (e^)=m + ^-^-((m + l)^-l)x fc . 



D a (J m (C))= [—) J m (e 2 



x=0 



□ 



Corollary 4.3. Suppose a, m are integers a > and m > 1. Then 

( 2(m a+1 - (m + l) a+1 )G a+1 (g) if a is odd, 
D a (WC,?))= lm+l i/a = 0, (4.22) 

I otherwise. 

Proof. The proof of (I4.22p when a = is straightforward. Suppose a is even and 
positive. Then by Lemma [4.21 

D a (L , m (C,q)) = D a (K 0tm (C)) 

a + 1 
= 0, 

since the Bernoulli numbers are zero when k > 1 is odd. Now suppose a is odd. 
Then again by Lemma 14.21 

A*(A),m(C, <?)) = £>a(#Q,m(C)) + 2(m a+1 - (m + l) a+1 )$ a (g) 

= ^±i(_ m «+i + ( m + + 2 (m a+1 - (m + l) a+1 )$ a (g) 

a + 1 

= 2(m" +1 -(m + l) a+1 )C7 Q+1 (g). 

□ 

By applying -D a -i to both sides of (14.151) and using ( 14. 22ft we obtain the following 
recurrence 

D a (F 0>m (C,q)) = (m + §) A*-! {F 0>m {C,q)) 

+ E 2 ( 2 ° _ 1 (m ' 4 ~ (m + l) 2i )G 2i (g)^(F , m (C, ?))■ (4.23) 
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This together with the initial value 
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A)(iV(C, q)) = F , m (l, q) = (4.24) 
uniquely determines the coefficients D a (F 0m ((,q)). We compute some examples 



£0(^0,2) 
£1(^0,2) 



3 

2' 
15 

T 

75 



£2(^0,2) = — - 15G 2 = 10 - 15$ 1; 
_ . 365 225^ 225 225^ 

L> 4 (F 0j2 ) = - ^-G 2 + A50G 2 2 - 195G 4 = 82 - 600$i + 450$? - 195$ 3 - 

■_)_ _ 



4.4. The terms Fj^ m ((, q) (1 < j < m — 1). Suppose 1 < j < m — 1. We may obtain 
a similar recurrence for D a (Fj t7n ((, q)). This time we find that 

S c Fj, m (C: q) = %m(C, q) Fj, m (C, q), (4-25) 

for some function Lj^ m ((,q) that satisfies 

{2 ((m + i + l) a+1 - (to - z) a+1 ) G a+1 (g) if a is odd, 
—j (to + |) if a = 0, (4.26) 

otherwise. 

The proof of (I4.26P is analogous to the proof of Corollary 14.31 By applying D a _i to 
both sides of (I4.25[) and using ( I4.26[) we obtain the following recurrence 



D a (F j>m ((,q)) = §)A*-i(*i,m(C,?)) 
L«/2J j / _ i \ 

+ E E 2 ( 2i- 1 ((m + * + 1)21 " (m - ^ 2i )^(g)^- 2i (Fi,m(C, <?))• 

(4.27) 



This together with the initial value 

A>(fi,m(C, ?)) = ^.m(l, g) = II J+Tll (428) 



m — i 



m + i + I 
i=i 
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uniquely determines the coefficients D a (Fj >m (£,q)). We compute some examples 
1 

~4' 
5 



D (F lt2 ) 



/ 7-i \ 25 15^ 5 15 

D 2 (F l2 ) = G 2 = $i, 

v ' ; 16 2 4 2 

„ /rn , 125 225^ 25 225^ 
D 3 (F 12 ) = 1 G 2 = 1 $i, 

D 4 ( Fl 2 ) = -— - — G 2 - 675C7 2 . - = - - 225$x - 675$? - — $ 3 . 

4v m; 64 4 2 2 4 1 2 

4.5. The terms 5 , 2 m +i(C"' +1 ) 2 ? ?) (1 < j < m — l). Again suppose that 1 < j < m—1. 
Consider the operator Tj that operates on a function f(() by Tj(f(()) = /(C J )- Then 

5 c o^=j (Tjo^). (4.29) 

A simple induction argument gives 

% o Tj = f (Tj o , (4.30) 

and 

D a oT 3 = 3 a D a . (4.31) 

Thus by (14. 10[) we have 

D e S 2m+1 (C +1 ,z,q) = (j + iyP 2m+1/ {W 2m+l )^ 2m+l \z,q). (4.32) 

4.6. The main theorem. We are now ready to derive our main theorem. Applying 
D 2m to both sides of flUD and using fl4TTU|) . iH, (Q7|l . (|Q5jl . and 
(gSp we have 

Theorem 4.4. 

(-l) m+1 (2m)! (m + 1)! (m - l)\[C*(z, q)? m+1 {qf™ +1 

P2m+l,2m(H 2m +l) 
j=l a =0 ^ ' 

- £> 2w (Fo,m(C, ?)) j S( 2m+1 ) (z, q) (4.33) 

where the coefficient functions D a (Fj^ m )((,q)) (0 < j < m — 1 ) are given recursively 
by h4-23\ ) and fl^.£7| ), and their initial values \\-2$ and h4-28\ ). 

For n > let V n be the Q-vector space spanned by the monomials with 
a + 2b + 3c = n. We Define 

n 

W n = ^V„; (4.34) 

fc=0 
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i.e., W n is the Q- vector space spanned by the monomials $"$^$5 with < a + 26 + 3c < 
n. We call W n the space of quasi-modular forms of weight less than or equal to 2n. 
This agrees with the definition in (2j p. 355] except this time we allow monomials of 
weight 0. 

Corollary 4.5. Suppose m > 1. Then there exist quasi-modular forms fj £ Wj for 
1 < j < m such that 

m— 1 \ 

Ci + E S(2m+1) (^?) = (2m)! [C*(,,g)] 2m+1 (g)L m+1 . (4.35) 

fc=0 / 

Proof. Suppose m > 1. It is well-known that 

$2n-l £ W„. 

See [2j Eq.(3.25)]. Equation (I4.19p . the recurrence (I4.27P and a simple induction argu- 
ment imply that 

D a (F j>m ((,q))eW la/2b 
for 1 < j ' < m — 1. Similarly using (I4.22p and (I4.23P we find that 

D2m{F , m {C,q))eW m . (4.36) 

Now we calculate the coefficient of T-L^m+i i n ^ ne light side of (I4.33p . The degree of 
the polynomial P2m+i,2m-a(%) is [(2m — a)/2j. Assuming < [(2m — a)/2j we have 
2A; < 2m — a and [_ a /2j < m — k, and in this case D a (Fj ym ((,q)) is in Wm-k- This 
together with (I4.36P implies that the coefficient of %2m+i ls i n W m -k for < k < m. 
The coefficient of "H^m+i i s 

/c = 2 + 2E (- 1 y c ,- + 1 )-n^ I 

j=l fc=l 

^s-to+d- , (m -;>;' m+1)! ,., (4.37) 

^ (m- j - l)\(m + j + 1)! 

' 2m 



2(m-l)!(m + l)! ^ / 
(2m)! >{3+ > V 



m — j — 1 



We show that 



/ = (-l) m+1 (m + l)!(m-l)!. (4.38) 
In view of (I4.37P this is equivalent to showing that 

m— 1 



(2m,. 

which we can rewrite as 

m— 1 



T E(-ir +i+1 (i+W 2m = 1, 



(4.39) 



by replacing j by m — j — 1 in the sum. 



2^(-l)i(m-jW . =(2m)!, (4.40) 

o— n \ J / 
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Since 

m— 1 /„ \ 2m /„ 

2m 



;— n \ J / a— ™ i i V 



j=0 - j=m+l 

it suffices to prove 



i=o 

Beginning with the elementary identity 

2m 



^ /2m \ 

^(-iy(m-j) 2m . =(2m)L (4.41) 

o— n V J / 



i=o V ^ / 

setting x = — and y = we obtain 

£ m j ( _ 1} i r -i = c -m (c _ 1)2m _ (4 42) 

i=o ^ J ' 

We apply /} 2 m to both sides of (14.421) and argue as in Section I4~2l to obtain (I4.4ip which 
completes the proof of (I4.38p . The final result (14.351) follows from (I4.33P by dividing 
both sides by fo and using ( 14.381) . □ 

4.7. Some examples. We illustrate Theorem 14.41 and Corollary 14.51 with some exam- 
ples. We show details of the calculations for the cases m — 1, 2. In cases m = 3, 4 we 
give the quasi-modular forms fj (1 < j ' < m) in Corollary 14.51 i n terms of the functions 
<&2k-i{q) rather than the G2k{q)- 
Example m = 1. 

4 [C*(z, q)} 3 ^ = ( 9 + 2 U 3 - D 2 (F 0tl (C q))\ (z, q) 
9 + 2U 3 - (5- 12$i)j£< 3 >(z, g ) 
2^3 + 4 + 12$! )rt 3) (z,q), 



and 



^3 + 2 + 6$! )^ 3 \z,q) = 2[C%z,q)} ;i (q) 



This identity implies the Rank-Crank PDE (11.71) as in [2], Section 2]. 
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Example m — 2. 

-lU[C*(z,q)f(q)t 

625 + 100 U\ + 2U* 2 5 + 16 D (F 1>2 ((, q)) (625 + 100 H\ + 2 U* 2 5 ) 
+ 32 D 1 (F 1)2 (C, g)) (125 + 15 H* 5 ) + 24 D 2 {F 1>2 {Q, q)) (25 + 2 H\) + 40 D 3 (F li2 (C, g)) 



+ 2 D 4 (F^(C, g)) - A(Fo, 2 (C, g))j s (5) (^, g) 

625 + 100 ?T 5 + 2 U*\ - 4 (625 + 100 H\ + 2 K*s) 

125 

+ 20 (125 + 15 7T 5 ) - (30 + 180 $0(25 + 2H\) + (— + 2250 $ x ) 
+ (i - 450 $i - 1350 - 255 $ 3 ) + (-82 + 600 $x - 450$? + 195 $ 3 ) j £ (5) (z, g), 
and 

jj + (60$x + 10) U\ + 300 $' + 10 $ 3 + 350 $ x + 24 J (2, g) = 24 [C*(z, g)] 5 (g)L- 

(4.43) 

In this case of Corollary 14.51 we see that 

h = 60$x + 10, 

f 2 = 300 $1 + 10 $ 3 + 350 $1 + 24. 
We show how this identity implies ()1.18p . We need the results 

5 g ((g)oc) = -$1 (g)oo, 6 q ($i) = |$! - 2$? + |$ 3 . 

This implies that 

Hl(V(z,q)) =Hl((q)lG®(z,q)) = (g)L(^ - 30$i)G< 5 ) (z, q), (4.44) 

and 

Hf(J:*(z,q))=Hl 2 ((q)lG^(z,q)) = (g)L(^5 2 -60$ 1 ^-50$ 1 + 1500$?-250$ 3 )G( 5 )(^ g). 

(4.45) 

Substituting (14T451) into (Q3l) we find 

{HI 2 + 10^ + 24 - 240$ 3 )G (5) (^ g) = 24 [C*{z, g)] 5 (g) 
which simplifies to (11.18P since 

n, = n; + 5. 



5 

00 > 



22 SONG HENG CHAN, ATUL DIXIT AND FRANK G. GARVAN 

Example m = 3. We find that 
h = 210 $i + 28 

f 2 = 210 $ 3 + 8820 + 252 + 4410 $i 

f 3 = 41160 + 2450 $ 3 + 14 $ 5 + 720 + 22736 $i + 2940 $ 3 $i + 102900 

Example m = 4. We find that 
f x = 504 $i + 60 

f 2 = 1260 $ 3 + 24948 $ a + 1308 + 68040 

/ 3 = 136080 $ 3 $i + 504 $ 5 + 45360 $ 3 + 2449440 + 403704 $ x + 2449440 + 12176 
/ 4 = 40320 + 2126232 $x + 404082 $ 3 + 9828 $ 5 + 2653560 $ 3 $ x + 21820428 

+ 9072 $x $ 5 + 47764080 + 18 $ 7 + 1224720 $ 3 + 11340 <$> 2 3 + 11022480 

5. Concluding remarks 

The main goal of this paper was to show how the generalized Lambert series identity 
(11.31) leads to the higher level Rank-Crank-type PDEs of Zwegers. The first author's 
proof |8j of (11.31) only involves a partial fraction argument and this together with the 
proof in Section H] gives an elementary g-series proof of these higher level Rank-Crank- 
type PDEs. The elliptic function proof of (11.31) in Section [2] is independent of the other 
sections. Our form of Zwegers's result (11.231) was given above in (I4.35p . In our form 
the coefficients are quasimodular forms rather than holomorphic modular forms. The 
quasimodular function E 2 occurs in Zwegers's result as part of the definition of his 
operator Tik- Our coefficient functions are given recursively. It would be interesting to 
find explicit expressions for the coefficients and to derive the form of Zwegers's result by 
our method. The coefficients in the 4 th order PDE (11. IS)) only involve the holomorphic 
modular form E4, and the differential operator H* does not involve the quasimodular 
E 2 . It would be interesting to determine whether there is a renormalization of higher 
order Rank-Crank-type PDEs which only involve holomorphic modular forms, either 
as coefficients or in the definition of the differential operator. Bringmann, Lovejoy 
and Osburn [5J, [6 J found Rank-Crank-type PDEs for overpartitions. Bringmann and 
Zwegers [7] showed how these results fit into the framework of non-holomorphic Jacobi 
forms and found an infinite family of these PDEs. However these PDEs only involve 
Appell functions of level 1 or 3. It would be interesting to determine whether the 
methods of this paper could be extended to find PDEs for higher level analogues. 
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